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Abstract. New results for the double beta decay of ^^Ge are presented. They are 
extracted from Data obtained with the Heidelberg- Moscow experiment, which oper- 
ates five enriched ^^Ge detectors in an extreme low- level environment in the Gran 
Sasso underground laboratory. The two neutrino accompanied double beta decay is 
evaluated for the first time for all five detectors with a statistical significance of 
47.7kgy resuhing in a half life of T^^i^ = [1.55±0.01(stat) lo:l5(syst)]x 10^^ years. The 
lower limit on the half-life of the 0/^/3/3 -decay obtained with pulse shape analysis is 
T'j'^a > l-9xl0^'' (3.1x10^^) years with 90% C.L. (68% C.L.) (with 35.5kgy). This re- 
sults in an upper limit of the effective Majorana neutrino mass of 0.35 eV (0.27 eV). No 
evidence for a Majoron emitting decay mode or for the neutrinoless mode is observed. 



1 Double Beta Decay 

There seems to be a general consensus over the neutrino oscillation interpretation 
of the atmospheric and solar neutrino data, delivering a strong indication for a 
non-vanishing neutrino mass. While such kind of experiments yields information 
on the difference of squared neutrino mass eigenvalues and on mixing angles, the 
absolute scale of the neutrino mass is still unknown. Information from double 
beta decay experiments is indispensible to solve these questions ||l|,|[||,^ . Another 
important problem is that of the fundamental character of the neutrino, whether 
it is a Dirac or a Majorana particle. Neutrinoless double beta decay could answer 
also this question. The HEIDELBERG-MOSCOW experiment is giving since 
almost eight years now, the most sensitive limit of all /3/3-experiments worldwide 
[||. Double beta decay, the rarest known nuclear decay process, can occur in 
different modes: 
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- decay : A{Z, N) ^ A{Z+2, N-2) + 2e- + 2De (1) 

- decay : A{Z, N) ^ A{Z +2, N -2) + 2e- (2) 
Oz^(2)x^/3 - decay : A{Z, N) ^ A{Z +2, N -2) + 2e- + {2)x (3) 

The two-neutrino decay mode (1) is a conventional second order weak process, 
allowed in the Standard Model of particle physics. So far it has been observed 
for about 10 different nuclei |^,^,D. An accurate measurement of the half-life of 
the decay is of importance, since it provides a cross-check on the reliability of 
matrix element calculations. The Majoron emitting decay mode (2) could reveal 
the existence of light or massless bosons, so called Majorons, with a non-zero 
coupling to neutrinos. The neutrinoless mode (3) is by far the most exciting 
one due to the violation of the lepton number by two units. It can not only 
probe a Majorana neutrino mass, but various new physics scenarios beyond the 
Standard Model, such as RjDarity violatingsupersymmetric models R-parity 
conserving SUSY models (7| , leptoquarks M , violation of Lorentz-invariance 
and compositeness (for a review see [|^ JlT| , |l2| , ^ ) . Any theory containing 
lepton number violating interactions can in principle lead to this process allow- 
ing to obtain information on the specific underlying theory. The experimental 
signature of the neutrinoless mode is a peak at the Q- value of the decay, whereas 
for the two-neutrino and Majoron-accompanied decay modes well defined con- 
tinuous energy spectra are expected. They are identified by their spectral index 
n, defined as the power of the energy in the phase space integral (see [Q). The 
Majoron emitting modes are characterized by n=l,3,7, while for the 2vj3(3 decay, 
n=5. 



2 The Heidelberg-Moscow experiment 

The Heidelberg-Moscow experiment operates five p-type HPGe detectors in the 
Gran Sasso underground laboratory which were originally grown from 19.2 kg of 
enriched ^^Ge. The total active mass of the detectors is 10.96 kg, corresponding 
to 125.5 mol of ^^Ge, the presently largest source strength of all double beta 
experiments. The enrichment of the used Germanium is 86%. A detailed de- 



scription of the experiment is given in |15 



To check the stability of the experiment, a couple of parameters such as 
temperature, nitrogen flow, leakage current of the detectors, overall and indi- 
vidual trigger rates are monitored daily. An energy calibration is done weekly 
with a ^^*Th and a ^^^Eu-^^*Th source. The energy resolution of the detectors 
at 2614 keV ranges from 3-3.7 keV. The energy thresholds for data recording are 
set to about 70keV (with exception of the second detector, which is used for 
dark matter measurements in addition, see [p^). 

Figure H shows the combined sum spectrum of all five enriched detectors 
of the Heidelberg-Moscow experiment with a statistical significance of 47.4 kg y 
(see |2^). The large peak-to-Compton ratio of the detectors facilitates the iden- 
tification of 7 activities. The easily identified background components consist 
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Fig. 1. Sum spectrum of all five Ge detectors after 47.4 kg y of measurement. The 
most prominent identified lines are labeled. 

of primordial activities of the natural decay chains from ^ssy g^^^^ 232 rpj^ from 
^°K, anthropogenic radio nuclides, like ^^^Cs, ^'^'^Cs, ^^^Sb and ^°^Bi and cos- 
mogenic isotopes, such as ^''Mn, ^^Co, ^*Co, ^"^Co and ^^Zn. Hidden in the 
continuous background are the contributions of the bremsstrahlungs spectrum 
of ^'^'^Bi (daughter of ^^°Pb), elastic and inelastic neutron scattering and direct 
muon induced events. 



3 Background model 

The evaluation of the spectra caused by the 2i//3/3 decay and the Majoron- 
emitting decay modes requires a detailed knowledge of the composition of the 
background on which they are superimposed. To unfold the background, a Monte 
Carlo simulation was performed. It is based on the CERN code Geant3.21, 
modified for simulating radioactive decays with the complete implemented decay 
schemes taken from . Five parts of the experimental setup have been identified 
to represent the main locations of the radioactive impurities: the LC2-Pb shield, 
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the copper shield, the copper and plastic parts of the cryostats and the Ge crys- 
tals themselves. Other materials or locations in the detector array are negligible 
due their small masses and low activities. The following background components 
were simulated: the natural decay chains of ^ssy ^j^j '^^'^Th, ^"K, cosmogenic 
and anthropogenic isotopes, muon showers and neutron induced interactions. It 
was assumed that the ^38^ g^j^j^ 232rpj^ decay chains are in secular equilibrium 
and that the radioactive isotopes in the respective materials are uniformly dis- 
tributed. Muon-induced showers were simulated based on the measured flux and 
energy distribution of muons in the Gran Sasso underground laboratory jl^ . Not 
considered were muon-induced neutrons in the detector shielding materials, due 
to still large uncertainties in the absolute n-flux determinations in Geant3.21. 
This component belongs to the non-identified background which will be discussed 
below. The measured neutron flux in the Gran Sasso underground laboratory 
[|9| was simulated using the MICAP implementation in Geant |Q . The activi- 
ties of and ^^"Pb in the LC2-Pb shield were determined in separate activity 
measurements In order to extract the best fit values for each activity, a 
least-squares method has been used. The location of the radioactive impurities 
was determined by comparing the peak intensities of multiline isotopes with the 
simulation. The error of a possible misplacement is part of the systematic error 
of the background model. The influence of each radioactive impurity located in 
one detector on all other detectors was considered. We identified a total number 
of 142 lines in the spectra of the five enriched Ge detectors. Their measured 
intensities were used to normalize the simulated components of the background 
model. Table ^ shows the identified background components, their estimated ac- 
tivities and their most probable locations in the experimental setup. The main 
background sources (natural decay chains, cosmogenics and anthropogenic ra- 
dionuclides) were located in the copper parts of the cryostats. In the Ge crystals 
themselves, only cosmogenic radionuclides were identified. There is no intrinsic 
U/Th contamination of the crystals, due to the absence of a-peaks in their high 
energy spectra (the single a- line at 5.3 MeV detected in two of the five detectors 
originates most likely from surface contaminations at the inner contact). Exter- 
nal a and /3 activities are shielded by the about 0.7 mm inactive zone of the 
p-type detectors on the outer crystal surface. Figure ^ shows the contribution of 
the simulated background components on the original measured sum spectrum 
of the Ge detectors (for details of the simulations see p^). 



4 Results for the 2v(3f3 and the 0^'(x)x/3/3 decays 

In Fig. ^the summed data of the five detectors are shown together with the result 
after subtracting the identified background components. A bin width of 20 keV is 
chosen in order to avoid statistical fiuctuations when subtracting the simulated 
7 lines from the measured spectrum. The contribution of the 2vf3(3 decay to 
the residual spectrum is clearly visible. Its half-life was determined under the 
assumption that the entire residual spectrum is composed of the 2vj3l3 -signal. 
Due to non-identified background in the energy region below 700 keV, the fit 
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Table 1. Identified background components (primordial, cosmogenic, anthropogenic), 
their estimated activities and most probable locations in the full setup of the 
Heidelberg-Moscow experiment. 



interval for the 2i//3/3-signal is chosen between 700-2040 keV. With the above 
assumption, this region contains 64553 2j//3/3 events, corresponding to 51.7% of 
the total 2vl3(3 -signal. 

The theoretically expected 2vj3(3 spectrum was fitted to the data in a 
maximum-likelihood-fit with as free parameter, resulting in the following 
half-life for the 2i//3/3 -decay at 68% C.L. (combined result for the five detectors): 

Tf = ( 1.55 ± O.Ol(stat) t^:}^(syst) ) x lO^^ y (4) 

The statistical error is evaluated from the parabolic behaviour of the log- 
arithmic likelihood ratio which corresponds to a function. The systematic 
error includes the error of the simulated detector response, the error made by 
the misplacement of background activities and the normalization error due to 
the statistical error of the measured 7-lines (see ||2^). 

The inferred value for the half-life is consistent with earlier results of this 
experiment |T5||2^ and with the result of j2^ , as well as with the range of the- 
oretical predictions, which lie between 1.5xl02°-2.99 xlO^^y ||2|,|6|,|^@ . 
The prediction of for the 2v[3l3 matrix element agrees within a factor of 

y/2 with the experimental value. 

The half-life limits of the Majoron-emitting decay-modes were determined 
from the same data set by fitting the 2vj3(3 and the QvxPP spectra simultane- 
ously. The considered Majoron models are described in |Q. Since the selected 
energy interval starts at 700 keV, an analysis of the decay-mode with the spectral 
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Fig. 2. The simulated background components (sliaded areas) compared witli tlie orig- 
inal measured sum spectrum for all five detectors. 
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Table 2. Half-life limits for the Majoron-emitting decay-modes and derived coupling 
constants using the matrix elements from jsij for different majoron models (n is the 
spectral index of the decay mode) . 



index n=7 (maximum at about 500 keV) was not possible. The results of the fits 
for n=l and n=3 are shown in Table ^ The 21/^(3 half-lifes extracted in the 
two-parameter fits are consistent within la with the exclusive double beta decay 
evaluation. In Table ^ a comparison of the effective Majoron neutrino couplings 
extracted for different double beta nuclei is made. 
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Fig. 3. Summed spectra of all five detectors after 47.7 kg y of measurement together 
with the residual spectrum after subtracting all identified background components. The 
thick line shows the fitted 2///3/3-signal. 
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Table 3. Half-life limits on the Majoron-emitting decay-mode GvxPP extracted from 
different nuclei and the derived limits on the effective Majoron-neutrino coupling for 
n=l. 



5 Results for the 0uf3(3 decay 

For the evaluation of the Qv[3(3 decay we consider the raw data of all five 
detectors as well as data with pulse shape analysis. The pulse shape analysis 
method used here is described elsewhere [p3|. No further data manipulation is 
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done, e.g. the previously established background model is not subtracted. We 
see in none of the two data sets an indication for a peak at the Q-value of 
2038.56±0.32keV Q of the Oiy/SP decay. 

The total spectrum of the five detectors with a statistical significance of 
53.9 kg y contains all the data with the exception of the first 200 d of mea- 
surement of each detector, because of possible interference with the cosmogenic 
^^Co. The interpolated energy resolution at the energy at the hypothetical 
peak is (4.23±0.14) keV. The expected background in the Oi/fSP region is es- 
timated from the energy interval 2000-2080 keV. In this range the background 
is (0.19±0.01) counts/ (kg ykeV). The expected background in the 3a peak in- 
terval, centered at 2038.56 keV interpolated from the adjacent energy regions, is 
(110.3±3.9) events. The number of measured events in the same peak region is 
112. To extract a half-life limit for the 0i^/3/3 -decay wc follow the conservative 
procedure recommended in 

With the achieved energy resolution, the number of excluded events in the 
3(T peak region is 19.8 (12) with 90% C.L. (68% C.L.), resulting in a half-life 
limit of (for the 0+ 0+ transition): 

> 1-3 X 10^5 y 90% C.L. 
T°^2 > 2.2 X 10^5 y 68% C.L. 

We consider now the data for which the pulse shape of each interaction of 
the detectors was recorded and analyzed. The total statistical significance is 
35.5 kg y and the background index in the energy region between 2000-2080 keV 
is (0.06±0.01) events/ (kg y keV) , about a factor 3 lower than for the full data 
set. This is due to the large fraction of multiple Compton scattered events in this 
energy region, which are effectively suppressed by the pulse shape discrimination 
method. The expected number of events from the background left and right of 
the peak region is (20.4±1.6) events, the measured number of events in the 3(t 
peak region is 21. Following again the method proposed by we can exclude 
9.3 (5.5) events with 90% C.L. (68 % C.L.). The limit on the half-life is: 

^1% > 1-9 X 10^5 y 90% C.L. 
T'i'/2 > 3-1 X 10^5 y 68% C.L. 

To examine the dependence of the half-life limit on the position of the 3a 
peak interval (12.7 keV) in the spectrum, wc shifted the peak interval between 
2028 keV and 2048 keV. It results in a variation of the half- life limit between 
2.5 X 10^^ y and 1.2 x 10^^ y at 90% C.L. (for the data with pulse-shape analysis). 
This demonstrates a rather smooth background in the considered energy region. 
Figure ^ shows the combined spectrum of the five detectors after 53.93 kg y and 
the spectrum of point-like interactions, corrected for the detection efficiency, 
after 35.5 kg y. The solid lines represent the exclusion limits for the two spectra 
at the 90% C.L. Using the matrix elements of |^ and neglecting right-handed 
currents, we can convert the lower half-life limit into an upper limit on the 
effective Majorana neutrino mass, which are listed in Table lil 
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energy [keV] 

Fig. 4. Sum spectrum of all five detectors with 53.9 kg y and SSE spectrum with 
35.5 kg y in the region of interest for the Ou/SP -decay. The curves correspond to the 
excluded signals with T^'j^ > 1-3 x W^^ y (90% C.L.) and Tf"^^ > 1.9 x 10^^ y (90% 
C.L.), respectively. 
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Table 4. Limits on the effective Majorana neutrino mass from the 0^(3/3 -decay of 
^^Ge calculated with the matrix elements from 1431. 



The HEIDELBERG-MOSCOW experiment is presently giving the most strin- 
gent upper limit on the Majorana neutrino mass, of 0.35 eV at 90% C.L. (0.27 eV 
at 68% C.L.). The values quoted in a previous paper jii), with a statistical sig- 
nificance of 24.2 kg y of data with pulse shape analysis, were 0.2 eV for the mass 
limit and 0.38 eV for the sensitivity of the experiment (both at 90% C.L.), after 
the recommendation of [Q. Thus, not unexpected, after additional 11.3 kgy of 
statistics, the limit on the effective neutrino mass approached the experimental 
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sensitivity, as defined in . The mass limit varies within a factor of less than 
two for different matrix element calculations (see the discussion in 0). 



6 Summary and Discussion 



We performed an analysis of the most recent data of the Heidelberg-Moscow dou- 
ble beta decay experiment. The data of the complete setup with five enriched 
''^Ge-detectors, with a total statistical significance of 47.4 kg y, were analyzed 
with respect to the two-neutrino and Majoron emitting decay modes for the 
first time. A Monte Carlo simulation based on a modified version of Geant3.21 
was performed in order to identify the most significant background sources and 
to establish a quantitative background model. The theoretical shapes of the 
2vj3[3 and QvxPP decay spectra were fitted in a maximum-likelihood fit to the 
resulting spectrum after subtraction of the background model from the mea- 
sured, summed spectrum of all detectors. The low-energy background of the 
HEIDELBERG-MOSCOW experiment requires further investigation. A possible 
background source not taken into account so far could be surface contaminations 
of the crystal and/or copper parts of the cryostats with ^^°Pb, which is produced 
and accumulated by the decay of ^^^Rn. This and other potential background 
sources will be implemented in a new Monte Carlo simulation using Geant4 
. A more complete background model will allow to determine the half- life of 
the 2vl3(3 decay with still higher precision. 
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lines) (see 0,|,|,@) 



The resulting half-life for the 2ul3[3 decay confirms our previous measure- 
ment and confirms theoretical expectations p5| , ^ within a factor of two (a 
factor of \/2 in the matrix element). No evidence for a Majoron-accompanied 
decay or for the neutrinoless decay was observed. The upper limit on the effec- 
tive Majorana neutrino mass of 0.35 eV (0.27 eV) (using the matrix elements of 
[^) is the worldwide most stringent limit up to now. In Fig. ^this value from 
the HEIDELBERG-MOSCOW experiment is compared with limits of the most 
sensitive other /3/3-experiments. With this result for the limit of the effective 
Majorana neutrino mass double beta experiments start to enter into the range 
to give a serious contribution to the neutrino mass matrix (Fig. H). 

In degenerate models we can conclude from the experimental bound an up- 
per limit on the mass of the heaviest neutrino. For the Large Mixing Angle 
(LMA) MSW solution of the solar neutrino problem we obtain r7ii,2,3 < 1.1 eV, 
implying J^i^i < 3.2 eV ||l|,3,H ■ This first number is sharper than what has re- 
cently been obtained from tritium (m < 2.2 eV), and the second is sharper than 
the limit X^i'™* < 5.5 eV still compatible with most recent fits of Cosmic Mi- 
crowave Background radiation and Large Scale Structure data (see e.g. Q). The 
present sensitivity of the HEIDELBERG-MOSCOW experiment probes cosmo- 
logical models including hot dark matter already now on a level of future satellite 
experiments MAP and PLANCK (see 0,^). It is of interest also for new Z-burst- 
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models recently discussed as explanation for super-high energy cosmic ray events 
beyond the GZK cutoff 

The result for < m > from the HEIDELBERG-MOSCOW experiment has 
found large resonance, and it has been shown that it excludes for example the 
Small Mixing Angle MSW solution of the solar neutrino problem in degener- 
ate scenarios, if neutrinos are considered as hot dark matter in the universe 
[ ^2[p^ , p^ . This conclusion has been drawn, before the Superkamiokande collab- 
oration presented their evidence for exclusion of SMA MSW solution, in June 
2000. 

If future searches will show, that rrii, > 0.1 eV, then the three-neutrino mass 
schemes, which will survive, are those with neutrino mass degeneracy, or four- 
neutrino schemes with inverse mass hierarchy (see Fig. ||and |p],p|j48| ) . A substan- 
tial increase in sensitivity of double beta experiments beyond this level, requires 
new experimental approaches, making use of much higher source strength and 
drastically reduced background. This could be accomplished by our proposed 
GENIUS project which, operating 0.1-10 tonnes of enriched ^^Ge directly in 
ultrapure liquid nitrogen, could test the effective Majorana neutrino mass down 
to 0.01 or even 0.002 eV. 
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